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A new framework (ADHYDRO) is presented that has been designed for description of 
highly-anisotropic systems produced at the very early stages of relativistic heavy-ion col- 
lisions. The structure of ADHYDRO is very much similar to perfect-fluid hydrodynamics, 
however, ADHYDRO includes dissipation effects. Dissipation is defined by the form of 
the entropy source that depends on the pressure anisotropy and vanishes for the isotropic 
systems. With a simple ansatz for the entropy source, that obeys general physical re- 
quirements, we obtain the system of equations describing isotropization of the anisotropic 
systems. 

1 Introduction 

The soft hadronic data obtained in the heavy- ion experiments at RHIC, in particular the mea- 
sured values of the elliptic flow coefficient i>2, are most often interpreted as the evidence for a 
very fast equilibration of the produced matter (most likely within a fraction of 1 fm/c) and for 
its perfect- fluid behavior [TJ [5] . Such very fast equilibration is naturally explained by the as- 
sumption that the produced matter is a strongly coupled quark-gluon plasma (sQGP) [31121 [5]. 
There are, however, other explanations that assume that the plasma is weakly interacting. In 
this case the plasma instabilities lead to the fast isotropization of matter, an effect that helps 
to achieve the full equilibration in a short time [BJ. 

Several recent calculations have demonstrated that the large values of the elliptic flow mea- 
sured at RHIC may be successfully reproduced in the models that do not assume fast equili- 
bration. In Ref. [7] the stage described by the perfect-fluid hydrodynamics was preceded by 
the free streaming of partons (see also [5JIS]), while in Refs. [TUJ HJ the authors assume that 
only transverse degrees of freedom are thermalized. The approach towards the full equilibration 
was discussed in this context in Refs. [T^[T3J[T3]. Such results indicate that the assumption 
of the fast equilibration/isotropization might be relaxed (see, however, the very recent results 
concerning the directed flow |15i 116] . which suggest that thermalization may be indeed very 
fast). 

We note that the idea of practically instantaneous equilibration seems to contradict the 
results of the microscopic models of heavy- ion collisions. These models typically use the concepts 
of color strings or color-flux tubes. The matter produced by strings is highly anisotropic; the 
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transverse pressure is usually much larger than the longitudinal pressure (the transverse and 
longitudinal directions are defined with respect to the beam). Similar situation takes place in 
the Color Glass Condensate (CGC) approach where the longitudinal momentum distribution 
is much narrower than the transverse one [T7J [T51 Q15] . 

In this talk we describe briefly the framework of ADHYDRO — highly-anisotropic and 
strongly-dissipative hydrodynamics [20j [21] . This framework has been designed for modeling of 
very early stages of heavy-ion collisions, where high anisotropy of the produced system is ex- 
pected. In Ref. [3D] ADHYDRO was used to describe boost- invariant, one-dimensional systems, 
whereas in Ref. [H] ADHYDRO was applied in non-boost-invariant situations. Sensitivity of 
the model to various modifications of its assumptions has been analyzed recently in [22] . 

The proposed model has a structure that is very much similar to the perfect-fluid hydrody- 
namics. The main two differences are connected with: i) the possibility that the longitudinal 
and transverse pressures arc different, and ii) the possibility of entropy production. We note 
that by relaxing the assumption about the isentropic flow, we have generalized our previous 
formulations of anisotropic (magneto) hydrodynamics presented earlier in Refs. [2"31 |2"4"] . 

It is important to emphasize that deviations from equilibrium are naturally described in the 
framework of viscous (Israel-Stewart) hydrodynamics, for a recent review see [25j . However, 
the region of the applicability of viscous hydrodynamics extends to the systems that are close 
to equilibrium (i.e., with small anisotropy). Thus, from the formal point of view, the use of 
viscous hydrodynamics in the description of very early stages of the heavy-ion collisions may 
be inadequate, since the strong reduction of the initial longitudinal pressure leads to significant 
deviations from equilibrium. 

One may show that for small pressure anisotropics and purely longitudinal expansion AD- 
HYDRO agrees with the 2nd order Israel-Stewart viscous hydrodynamics (we note that longi- 
tudinal expansion dominates at the early stages of heavy- ion collisions). On the other hand, for 
very large anisotropies ADHYDRO should be treated as a hydrodynamics-like model, that has 
no direct relation to the Israel-Stewart theory. Certainly, if the anisotropies are very large, one 
should use the kinetic theory (which has its own well known limitations) or consider specific 
models of the collisions. ADHYDRO is a proposal for such modeling. Other possibility is to use 
the method presented in Refs. [15] [16] [26] , where an exponentially decreasing time dependence 
of the difference P± — Pii is assumed. 

A work along similar lines have been recently done by Martinez and Strickland [27] |28| , who 
derive analogous equations from the Boltzmann equation with the collision term treated in the 
relaxation time approximation. Instead of the usual expansion around the isotropic one-particle 
distribution, they consider an expansion around an anisotropic distribution. 

2 Basic equations 

The ADHYDRO model is defined by the following form of the energy-momentum tensor (23] [24] 

= (e + P±)U> t U u -P±gT ~ (P± ~ P\\)V*V U , (1) 

where e, P_l, and Pi are the energy density, transverse pressure, and longitudinal pressure, 
respectively. In the case Pj_ = P\\ = P the form of the perfect-fluid energy-momentum tensor 
is reproduced. The four-vector [/ M describes the standard hydrodynamic flow 

U»=y(l,v x ,v y ,v z ), 7 =(l-« 2 r 1/2 , (2) 
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and V 1 - 1 defines the longitudinal direction (corresponding to the collision axis) 



V" = 7,(^,0,0,1), lz = {l-v 2 z ) 



-1/2 



(3) 



The four-vectors t/ M and V 1 satisfy the following normalization conditions: U 2 = 1,V 2 = —1, 
and U ■ V = 0. In the local-rest-frame of the fluid element, where v = 0, the four-vectors U 11 
and are reduced to the simple forms, = (1,0,0,0), = (0,0,0,1), and the energy- 
momentum tensor has the following diagonal structure 
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In addition to the energy-momentum tensor ((T|) we introduce the entropy flux 

o-" = all", 



(4) 



(5) 



where a is the entropy density. The quantities 
and P\\ . In particular, for massless partons 



and a are functions of the two pressures, P± 
II . in particular, lor iiiassiess partoiis i r ^ ~ and e = 2P± + P|| . 
The space-time dynamics of the system is described by the following equations 



= 0, 
<V = E. 



(6) 
(7) 



Equation ([6]) expresses the energy-momentum conservation, while Eq. ([7]) describes the entropy 
production. The form of the entropy source E determines the dynamics of the anisotropic 
system. In Ref. [20] we proposed the following ansatz 



E = 



(i-v^) 2 



(8) 



where x is the parameter describing the anisotropy of the momentum distribution function of 
partons and r oq is a time-scale parameter. One may check that x ~ (P±/ P||) 4 / 3 [20] . 

The expression on the right-hand-side of Eq. © has several appealing features: i) it is 
positive, as required by the second law of thermodynamics, ii) it has a correct dimension, iii) it 
vanishes in equilibrium, where x = 1. Moreover, for small deviations from equilibrium, where 
\x — 1| <C 1, we find 

^ (*-!) 2 _ 



4t, 



(9) 



cq 



The quadratic dependence displayed in ^ is characteristic for the 2nd order viscous hydrody- 
namics where the entropy production is proportional to the viscous stress squared, E ex II 2 [29] . 
As is shown in |27j . II oc x — 1 in our case. 



3 Results for boost-invariant systems 

If the functions a(P± , P\\ ) and E(Pj_ , P|| ) are specified, Eqs. (|6]) and ([7]) form a closed system of 
5 equations for 5 unknown functions: three components of the fluid velocity v, the transverse 
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pressure P± , and the longitudinal pressure Pn . In the case of the longitudinal and boost- 
invariant evolution only two functions are independent. Following Ref. }20) . we express all 
thermodynamics-like quantities in terms of x and a rather than in terms of P± and Py . The 
four-vectors and are defined by the expressions 

U» = (cosh?7,0,0,sinhr?), (10) 
V 1 = (sinh?j, 0,0, cosh r?), (11) 

where r\ is the space-time rapidity, 77 = l/21n(t + z)/(t — z). Substituting Eqs. (fit)]) and (fTTj) 
in Eqs. © and ([7]) we obtain 

de(a, x) _ e(a, x) + P\\ (a, x) 



dr 
da a 

dr t 



V{v,x), (13) 



where r = \/t 2 — z 2 is the proper time. The explicit forms of the functions e{a, x) and P11 (a, x), 
determined by the microscopic phase-space distribution function, have been given in 20 . 

We solve Eqs. (fT2"j) and (JTSJ) numerically with x = xo set at t = To = 0.2 fm. The results of 
the microscopic models suggest that P|| <?C P± at the early stages of the collisions, hence we con- 
sider the case xo = 100. The time evolution is studied in the time interval: 0.2 fm < r < 10 fm. 
In Fig. [1] we show the time dependence of various physical quantities obtained for three differ- 
ent choices of the relaxation time: r oq = 0.25 fm (solid line), r cq = 0.5 fm (dashed line), and 
T cq = 1.0 fm (dotted line). Figure Q] (a) shows the time dependence of the asymmetry param- 
eter x. We observe fast changes of x at the beginning of the evolution, which are followed by 
more moderate changes at later times. The behavior shown in Fig. [T] (a) indicates that x ~ 1 
for t > 2r Gq . This is a desired effect showing that the system approaches the state of local 
equilibrium. 

In Fig. Q] (b) we compare the time evolution of the entropy density obtained from Eq. (|13j) 
with the Bjorken solution obj = uoTq/t. Here cto is the initial value of the entropy density. 
We note that the specific value of <7o is irrelevant for our analysis. The amount of the entropy 
produced in the regime described by the anisotropic hydrodynamics depends in our case on 
the relaxation time. For r oq = 0.25, 0.5, and 1.0 fm the entropy increases by about 70%, 85% 
and 105%, respectively. For r 3> r oq the ratio (ctt) / '{<JqTo) saturates indicating that the flow 
attains the form of the Bjorken flow. This behavior shows explicitly that our framework may 
be used to model the transition between the highly anisotropic initial behavior and the perfect- 
fluid stage. In Fig. [T] (c) we show the time dependence of the ratio of the longitudinal and 
transverse pressure. Again, we show three different time evolutions corresponding to three 
different relaxation times. For r ^> r oq the ratio approaches unity and the two pressures 
become equal. 



4 Summary 

We have introduced the new framework of highly-anisotropic hydrodynamics with strong dissi- 
pation. The effects of the dissipation are introduced by the special form of the internal entropy 
source. The source depends on the pressure anisotropy and vanishes for the isotropic systems 
to guarantee that the perfect-fluid behavior is reproduced for the locally equilibrated system. 
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Our numerical results have been presented for a simple one-dimensional system. Neverthe- 
less, the proposed formalism is general and may be applied to more complicated 2+1 and 3+1 
situations. In addition, different forms of the entropy source inspired by different microscopic 
mechanisms may be analyzed. 

The results presented in this paper were obtained in collaboration with Radoslaw Ryblcwski. 
This work was supported in part by the Polish Ministry of Science and Higher Education, grant 
No. N N202 263438. 
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Figure 1: (a) Time dependence of the asymmetry parameter x for three different choices of 
the relaxation time: r oq = 0.25 fm (solid line), r cq = 0.5 fm (dashed line), and r cq = 1.0 fm 
(dotted line), (b) Entropy density divided by the corresponding values obtained in the Bjorken 
model, (c) Ratio of the longitudinal and transverse pressures shown as a function of the proper 
time. All results are obtained with the initial asymmetry xq = 100. 
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